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Abstract - The Geosynchronous Imaging Fourier Transform 
Spectrometer (GIFTS) Sensor Module (SM) Engineering 
Demonstration Unit (EDU) is a high resolution spectral imager 
designed to measure infrared (IR) radiance using a Fourier 
transform spectrometer (FTS). The GIFTS instrument 
gathers measurements across the long-wave IR (LWIR), 
short/mid-wave IR (SMWIR), and visible spectral bands. The 
raw interferogram measurements are radiometrically and 
spectrally calibrated to produce radiance spectra, which are 
further processed to obtain atmospheric profiles via retrieval 
algorithms. This paper describes the processing algorithms 
involved in the calibration. The calibration procedures can be 
subdivided into three categories: the pre-calibration stage, the 
calibration stage, and finally, the post-calibration stage. 
Detailed derivations for each stage are presented in this paper. 
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1. INTRODUCTION 


2. CALIBRATION ALGORITHMS 
2.1 Pre-calibration Procedures 

In the pre-calibration stage, a phase correction algorithm is applied 
to the decimated and filtered complex interferogram I m (n). Its 
complex spectrum N m (o) has the fonn of 


N m {(T) = N((T)e J ^ a) +n(a), (1) 


where n(a) is the noise vector, and N(a) is the real and noise-free 
spectrum. To estimate and correct the phase function </>(o), we 
adopt the Forman-Vanesse-Steel (FVS) method (Forman, 1966). 
In this algorithm, a Hamming window w(n) is first applied to the 
complex interferogram with respect to the zero path difference 
(ZPD) at n—n zpc i, i.e., 


G(») = Hn)I m (n), (2) 


The Geosynchronous Imaging Fourier Transfonn Spectrometer 
(GIFTS) Sensor Module (SM) Engineering Demonstration Unit 
(EDU) is a high resolution spectral imager designed to measure 
infrared (IR) radiance using a Fourier transform spectrometer 
(FTS). The GIFTS instrument employs three Focal Plane Arrays 
(FPAs), which gather measurements across the long-wave (685- 
1130 cm' 1 ) IR (LWIR), short/mid-wave (1650-2250 cm' 1 ) IR 
(SMWIR), and visible spectral bands (Best, 2001). The raw 
interferogram measurements are radiometrically and spectrally 
calibrated to produce radiance spectra, which are further processed 
to obtain atmospheric profiles via retrieval algorithms. 

The calibration procedures can be subdivided into three categories. 
In the pre-calibration stage, a phase correction algorithm is applied 
to the decimated and filtered complex interferogram. The 
resulting imaginary part of the spectrum contains only the noise 
component of the uncorrected spectrum. The phase correction 
operation is perfonned on a set of interferogram scans for the 
ambient, hot blackbody references, and scene measurements. In 
the calibration stage, we first compute the spectral responsivity 
based on the previous results, from which, the calibrated ambient 
blackbody (ABB), hot blackbody (HBB), and scene spectra can be 
obtained. In the post-calibration stage, we estimate the Noise 
Equivalent Spectral Radiance (NESR) from the calibrated ABB 
and HBB spectra. We then implement a correction scheme that 
compensates for the effects of fore-optics. The FPA off-axis effect 
correction technique also is presented. 
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where W is the odd-numbered window length. The windowed 
interferogram l„/n) is circular shifted to obtained I’ w (n) such that 
/ \/0) =I w (n zpd ). The estimated phase function can be derived from 
the spectral phase response of the windowed interferogram, which 
is given by 


N w {cj) = F\[' w (n)\ 
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1 Im(A\,,(o-)) 
Re(V w (o-)) 


(4) 

( 5 ) 


where F{.} denotes the Fourier transfonn operation. Next, we 
apply the estimated phase function to the complex spectrum: 


Corresponding author. 



N(a) = N m (a)e 

= N(a)e^ a) -^ a) Kn(a)e~J^ G \ 

The resulting real and imaginary parts of the spectrum are 

Re{N(<j)} = N(<r) + Re{n(<j)}cos(0(<T)) 

Im{A(<y)} = Im{n(<j)}sin(^((j)). 

Note that the imaginary part contains only the noise component of 
the uncorrected spectrum in Equation (7). 


Figure 1(b) illustrates the phase responses for HBB and ABB; here 
the mean phase responses within the spectral band of interest are 
approximately zero. 

The phase correction operation is performed on a set of 
interferogram scans for the scene measurements and internal 
blackbody calibration references at ambient and hot temperatures. 
The resulting real spectra can be written as Nj H (o), N, A (o '), and 
N/ S (g) for HBB, ABB, and scene, respectively, where i denotes the 
scan number. 

2.2 Calibration Procedures 

In the calibration stage, the spectral responsivity, based on the 
previous results and ideal Planck blackbody spectra at given 
temperatures, can be computed from 
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in which, 
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represent the mean HBB and ABB spectra over S scans, 
respectively. Using the relation of 
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Figure 1. FW spectral responses after the phase correction for 
Pixel (64, 65): (a) real and imaginary spectra for HBB (286K) and 
ABB (260K), and (b) phase responses for HBB and ABB. 

Figure 1 depicts the spectral response functions after the phase 
correction procedure for a FW pixel element over 25 scans. In 
Figure 1(a), the real and imaginary spectra are shown. Note that 
the imaginary spectra contain only noise after the correction. 



N H (ct) = R(g)(B H (a) + 0(a)), 
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The offset term can be solved by 
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The calibrated HBB, ABB, and scene spectra, therefore, can be 
determined from 
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respectively, for i=l,2,...,S. Figure 2 illustrates the LW HBB, 
ABB, and scene spectra obtained from the calibration procedure. 



Since the calibration results described in Section 2.2 were 
obtained based on internal blackbody references, it is necessary to 
implement a correction scheme that compensates for the effects of 
fore-optics. This task can be achieved by using additional data 
measurements collected from an extended source. Given 

IT A 

Cy(cr)and C g (cr) as the calibrated extended blackbody sources 

averaged over all scans at hot and ambient temperatures, 
respectively, we have 

C e H (a) = R e (a)Bj/(c7) + O e (c7), (17) 

C e A (a) = R e (a)B^(a) + O e (a). (18) 

By solving RJ a) and O e (a), the actual scene radiances can be 
approximated as 


Figure 2. The calibrated LW HBB, ABB, and scene spectra 
for Pixel (64, 65) over 25 scans. 

2.3 Post-calibration Procedures 

In the post-calibration stage, the Noise Equivalent Spectral 
Radiance (NESR) can be estimated from the calibrated ABB and 
HBB spectra. The NESR is generally considered as a measure of 
the instrument noise performance, and can be estimated as the 
standard deviation of calibrated radiance spectra from multiple 
scans. The definitions for NESRs are given in Equations (15) and 
(16). Figure 3 shows the smoothed HBB and ABB NESR 
estimates over 25 scans for a single pixel. 
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Figure 4 shows the calibrated LW radiance spectrum before and 
after performing the fore-optic effect correction procedure. 
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Figure 4. The calibrated LW scene spectrum before and 
after the fore -optic effect correction. 

The final step in the post-calibration process involves the 
correction of the FPA off-axis effect. Since the GIFTS FPAs 
contain 128x128 detector elements, the distortion caused by the 
off-axis effect is too great to be neglected. The off-axis pixels are 
sampled at slightly shorter optical path differences (OPDs), which 
cause the spectra of these off-axis pixels to expand to higher 
wavenumbers (Kauppinen, 2001). The off-axis correction can be 
fonnulated as a fractional sampling rate conversion rate problem, 
which can be solved via sine interpolations, i.e., 


Figure 3. The smoothed LW NESR estimates for HBB and 
ABB (Pixel (64, 65), 25 scans). 
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where a is the on-axis wavenumber scale, a’ is the off-axis scale, 
and / is the off-axis factor that can be obtained from the FPA 
geometry. However, the sine interpolation method is 
computational intensive, furthermore, the truncation window 
associated with the sine kernels may cause overshoot in the 
resulting data. It has been shown that the zero-padding 
interpolation in the interferogram domain is equivalent to the sine 
interpolation in the spectral domain except that the zero-padding 
method cannot be evaluated at an arbitrary point (La Riviere, 
1998). To overcome this limitation, we designed the method of 
“over-padding”, in which, we assign an over-padding factor g that 
is closely associated with the final resolution of the corrected 
spectrum. By over-padding the interferogram g times, a fractional 
sample rate conversion can be realized. For instance, if f=0.9977, 
and the length of the spectrum is N, then the zero-padded 
interferogram length is the round off of g*N/f If g=100, the 
actual off-axis factor of 0.997702 is obtained; if g=l, then the 
actual value of /becomes 0.99758 due to round off errors. The 
corrected spectrum can be obtained by downsampling the Fourier 
transform of the over-padded interferogram by the factor g. 
Figure 5(a) depicts the calibrated radiance spectra for two LW 
pixels before any off-axis effect correction was perfonned. One 
pixel is located near the center of the FPA while the other is 
situated near the corner of the FPA. The wavenumber shift can be 
seen between the two spectra, especially for higher wavenumbers. 
In Figure 5(b), the corrected spectra for these two pixels are 
shown. An over-padding factor of g=20 was assigned during the 
correction. As a result, a more accurate alignment was achieved 
between these two spectra. 




(b) 

Figure 5. The off-axis effect correction: (a) center and comer 
pixels before the off-axis correction, and (b) center and comer 
pixels after the off-axis correction with the over-padding factor set 

to 20. 

3. CONCLUSIONS 

This paper describes the three stages for calibrating the 
GIFTS SM EDU data spectrally and radiometrically. 

The first stage involves correcting the complex interferogram for 
its phase errors via the FVS algorithm. In the second stage, the 
spectral responsivity is computed, followed by the derivation of 
calibrated radiance spectra for HBB, ABB, and scene. In the last 
stage, the NESRs are estimated based on the results from the 
previous stages. We then compensate for the fore-optic effects 
using an extended source as the reference. A correction method 
for the off-axis effect is also presented. 
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